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oj ; Abstract 



Classical transport theory is employed to analyze the hot quark-gluon plasma 
at the leading order in the coupling constant. A condition on the (covariantly 
conserved) color current is obtained. From this condition, the generating 
functional of hard thermal loops with an arbitrary number of soft external 
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bosonic legs can be derived. Our approach, besides being more direct than 

On 

r~| ■ alternative ones, shows that hard thermal loops are essentially classical. 
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It is the purpose of the present Letter to derive the hard thermal loops (HTL's) of QCD 
from classical transport theory. The generating functional of HTL's (with an arbitrary 
number of soft external bosonic legs) arises as a leading-order effect in this formalism. 

In the diagrammatic approach to high-temperature QCD it was realized that a resum- 
mation procedure is necessary in order to take into account consistently all contributions at 
leading-order in the coupling constant |l| . Such contributions arise from one-loop diagrams 
with external soft legs and hard internal momenta, the so-called "hard thermal loops" . The 
HTL approach provides gauge invariant results for physical quantities. An effective action for 
HTL's was given in || by exploring the gauge invariance condition for the generating func- 
tional. This condition has been identified with the equations of motion for the topological 
Chern-Simons theory at zero temperature 0. Furthermore, the eikonal for a Chern-Simons 
theory has been used to obtain a non-Abelian generalization of the Kubo formula [|J , which 
governs, through the current induced by HTL's, the response of a hot quark-gluon plasma. 

Another successful description of the hard thermal loops in QCD is based on (a truncation 
of) the hierarchy of Schwinger- Dyson equations 0; the generating functional for HTL's is 
obtained by performing a consistent expansion in the coupling constant. Alternatively, one 
can start from the composite effective action by requiring its stationarity J7J, and by 
using the approximation scheme developed in ||. The gauge invariance condition for the 
generating functional can then be obtained as a consequence of symmetry properties due to 
the form of the momentum integration measure (see [0]). 

The diagrammatic techniques of JT|, as well as the consistent expansion in the coupling 
constant presented in ]5], although remarkably insightful, are technically very involved and 
necessitate lengthy computations. Furthermore, they are puzzling in respect of the very 
nature of hard thermal loops. One wonders to what extent these are quantum effects: on 
the one hand HTL's emerge in loop diagrams and are described by the Schwinger- Dyson 
equations of quantum field theory. On the other hand it is generally believed that hard 
thermal effects are UV-finite since they are due exclusively to thermal scattering processes. 

This motivates us to find an alternative, classical formalism for hard thermal loops in 
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QCD. The obvious setting for such a search is the classical transport theory of plasmas 
(see for instance ||). Our effort was encouraged by the fact that for an Abelian plasma of 
electrons and ions, the dielectric tensor that can be computed from classical transport 
theory is the same as the one that can be extracted from the hard thermal corrections to 
the vacuum polarization tensor Moreover, the same situation is encountered for non- 

Abelian plasmas [[K],|Il[]. The transport theory for colored plasmas has been established in 
|i0| . There has not been, to the best of our knowledge, any attempt at deriving the complete 
set of HTL's from this formalism. It is our purpose to do so, for the HTL's with an arbitrary 
number of external soft bosonic legs. 

Consider a particle bearing a non-Abelian SU(N) color charge Q a , a = 1,...,N 2 — 1, 
traversing a worldline x a (r), where r denotes the proper time. The dynamical effects of the 



spin of the particles shall be ignored, as they are typically small. The Wong equations (12 
describe the dynamical evolution of x M , p^ and Q a : 

m — — = £r , (la) 
dr 

m d £ = gQ a F^ Pv , (lb) 
dO a 

m^- = -gf abc p»A b l Q c . (lc) 

The f abc are the structure constants of the group, F£ v denotes the field strength, g is the 
coupling constant, and we set c = h = ks = 1 henceforth. Equation (jlbf ) is the non-Abelian 
generalization of the Lorentz force law, and eq. fllcl ) describes the precession in color space 
of the charge in an external color field A^. It is noteworthy that the color charge is itself 
subject to dynamical evolution, a feature which distinguishes the non-Abelian theory from 
electromagnetism. 

The usual (x, p) phase space is enlarged by including color degrees of freedom for colored 
particles. Physical constraints are enforced by the choice of the volume element dxdPdQ. 
The momentum measure 

d P = -^-26(p Q ) 5{p 2 -m 2 ) (2) 



(27T) 



guarantees on-shell evolution and positivity of the energy. The color charge measure enforces 
the conservation of the group invariants, e.g., for SU(3), 

dQ = d 8 Q 5(Q a Q a - q 2 ) 5(d abc Q a Q b Q c - q 3 ) , (3) 

where the constants q 2 and q% fix the values of the Casimirs. 

The one-particle distribution function f(x,p,Q) is the classical probability density for 
finding the particle in the state {x,p,Q}. This probability density evolves in time via a 
transport equation, 

m df{x £ Q) =C[f](x,p,Q), (4) 

where C[f](x,p, Q) denotes the collision integral, which we henceforth set to zero. Using 
the equations of motion ([I]), in the collisionless case, (f|) becomes the Boltzmann equation, 

f(x,p,Q) = 0. (5) 

A self-consistent set of (non-Abelian Vlasov) equations for the distribution function and 
the mean color field is obtained by augmenting the Boltzmann equation with the Yang-Mills 
equations: 

[D v F vtM ] a (x) = r a (x) . (6) 

The covariant derivative is defined as D^f = d^5 ac + g f abc A b . The total color current J fia (x) 
is given by the sum of all contributions from helicities and particle species, 

J^)= E E 3" a {x) = 2 J2 f a {*)- (7) 

helicities species species 

Each j fia (x) (spin and species indices are implicit) is computed from the corresponding 
distribution function: 

j»°( x ) = g J dPdQ p»Q a f(x,p,Q) (8) 
and it is covariantly conserved 
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(D,f) a (x) = , (9) 

as can be checked by using the Boltzmann equation. For later convenience, we define the 
total and individual current momentum-densities: 

J"°(x,p) = E E f a &P) , f a (x,P) = 9 JdQ p»Q a f(x,p,Q) . (10) 

helicities species 

The above framework is now employed to study the soft excitations in a hot, color- 
neutral quark-gluon plasma. In the high-temperature limit, the masses of the particles can 
be neglected and shall henceforth be taken to vanish. The wavelenght of the soft excitations 
is of order 1/g \A\, and the coupling constant g is assumed to be small. We then expand the 
distribution function f(x,p,Q) in powers of g: 

f = f i0) +9f {1) +9 2 f {2) + ..., (11) 

where is the equilibrium distribution function in the absence of a net color field, and is 
given by: 

f {0 \po)=Cn B , F (p ) . (12) 

Here C is a normalization constant and ub,f(po) = l/(e^' po ' =F 1) is the bosonic, resp. 
fermionic, probability distribution. 

At leading-order in g, the induced current (|10"D is 

j" a (x,p) = g 2 J dQp»Q a fV(x,p } Q) , (13) 

while the Boltzmann equation ([5]) reduces to 

P" {^~Sf abC KQ^ f {1) (x,p,Q)=p^Q a F^f^( Po ) . (14) 

Due to the softness of the excitation, the d x term in the above equation is of order g\A\, so 
we are taking into account consistently all contributions in order g. 

The equations ([□]) and ( |T4D yield the following constraint on the induced current: 
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[p-Dftx, P )r = g 2 ^p v Fl p ^(jdQQ a Q b f^(p )) , (15) 

where, from color symmetry, we expect that / dQ Q a Qbf <y0 \po) = Cb,f tib,f{Po) $1 w hh 
Cb = N, Cp — | for gluons, resp. fermions. Thus, upon summation over all species (N F 
quarks, Np antiquarks and one [(N 2 — l)-plet] gluon) and helicities (2 for quarks-antiquarks 
and massless gluon), (|T5| ) yields, 



[p-D J^x,p)] a = 2g 2 p^F^(x)^-{Nn B (p ) + N F n F (p )] . (16) 

apo 

This equation is equivalent^] to the condition (3.13) of ]F| for the induced current (this 
condition has been previously obtained in ||). 

The subsequent steps which lead to the gauge invariance condition for HTL's can be found 
in |7]]. Assuming that the induced current can be expressed as the functional derivative of 
an effective action allows one to derive, from (|TBp, a condition on the latter action. The 
result can be presented as (eq. (3.26) in jTJ): 
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4V2n 3 m 2 D d A + , (17) 



where d+ = C&d„ A + = Q$_A„ with Q% = ^(1, £) , and m D =gT^/^Ell is the 
Debye screening mass. The effective action generating the HTL's has the form (first derived 
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^ / dkA&x)A&x) - J ^yW{A + ) , (18) 

where dfl denotes integration over all angular directions of the unit vector ^. The gauge 
invariance condition ([[7]) has been solved for W 10,0], thereby obtaining closed expressions 
for the generating functional Thtl- This concludes our derivation of hard thermal loops 
from classical transport theory. 



1 Owc notational conventions differ from those of the e)(p 2 )-function appearing in (3.13) there 
is here contained in the (massless limit of the) momentum measure dP @. 
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Let us now summarize and discuss our results. The classical transport theory for a 
hot collisionless plasma of colored particles is a well-suited formalism for studying the hot 
quark-gluon plasma. The recently proposed HTL's arise naturally in this formalism. The 
spin-statistics theorem is the only "quantum mechanical" input we used in our analysis. Its 
use is entirely justified by the conditions of temperature and density which prevail in a hot 
quark-gluon plasma. With respect to alternative descriptions our approach is much simpler. 
It also demonstrates explicitly that the HTL's of QCD are a classical effect. 
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